Logarithmically modified scaling of temperature structure functions in thermal 
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Using experimental data on thermal convection, obtained at a Rayleigh number of 1.5 XlO , it 
is shown that the temperature structure functions (ATp), where AT r is the absolute value of the 
temperature increment over a distance r, can be well represented in an intermediate range of scales 
by r (,p ip(r) p , where the £ p are the scaling exponents appropriate to the passive scalar problem in 
hydrodynamic turbulence and the function <p(r) = 1 — a(lnr /ru) 2 ■ Measurements are made in the 
midplane of the apparatus near the sidewall, but outside the boundary layer. 
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A deep similarity between statistical properties of tur- 
bulent passive scalar fluctuations and those of tem- 
perature fluctuations in the turbulent thermal convection 
|2j has been discovered recently from numerical simula- 
tions (see, for instance, refs. 0,0 and the literature cited 
therein). In the present note, we will use experimental 
data and a brief theoretical reasoning to go further in 
this same direction. 

We consider turbulent convection in a confined con- 
tainer of circular crosssection and 50 cm diameter. The 
aspect ratio (diameter/height) is unity. The sidewalls are 
insulated and the bottom wall is maintained at a constant 
temperature, which is slightly higher than that of the top 
wall. The working fluid is cryogenic helium gas. By con- 
trolling the temperature difference between the bottom 
and top walls, as well as the thermodynamic operating 
point on the phase plane of the gas, the Rayleigh number 
(Ra) of the flow could be varied between 10 7 and 10 15 . 
We measure the temperature fluctuations T(t) at various 
Rayleigh numbers towards the upper end of this range, 
in which the convective motion is turbulent. The results 
used in this paper correspond mainly to Ra — 1.5 x 10 11 . 
Time traces of fluctuations are obtained at a distance of 
4.4 cm from the sidewall on the center plane of the ap- 
paratus. This position is outside of the boundary layer 
region for the Rayleigh numbers considered here. More 
details of the experimental conditions and measurement 
procedure can be found in ref. 0. 

For the structure functions (AT?), where the absolute 
values of temperature increment AT\T(t + r) — T(t)\ is 
evaluated as a function of the time interval r and the 
angle brackets indicate time averaging, we do not observe 
any clear scaling (see fig. 1). However, the normalized 
structure functions show unambiguous scaling as 
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FIG. 1: Temperature structure functions (AT*) against r 
(in log-log coordinates). Here and for other figures, Ra — 
1.5 x 10". 



It should be noted that temperature structure func- 
tions in thermal convection were investigated in ref. 
using the experimental data described in 0. In par- 
ticular, the data used in were obtained at the cen- 
ter of the convection cell. While the structure func- 
tions (AT^) themselves showed no scaling properties, two 
scaling ranges were observed for the normalized struc- 
ture functions S p (t). In the present paper, we study 
the data obtained near sidewall (see above), where suffi- 
ciently strong convection wind is present. For this case, 
we observe only one scaling interval for S p (t). 

Since, for a passive scalar 9 advected and diffused by 
hydrodynamic turbulence, the scaling of structure func- 
tions follows (in both experiments and direct numerical 
simulations) the relation 



(see fig. 2). The scaling exponents z p , extracted as slopes 
of the best straight-line fits shown in fig. 2 (the scaling 
range is from 0.1s to 10s), indicate multiscaling. These 
exponents are shown in fig. 3 as circles. 
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the corresponding exponents of the normalized structure 
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FIG. 2: Normalized temperature structure functions 
(AT?)/{AT?) p/2 against r (in log-log coordinates). Straight 
lines (the best fits) are drawn to indicate the multiscaling (2). 

functions can be readily calculated from the definition as 

z P = Cp - |Ca- (4) 

We also show in fig. 3 the values of z p , calculated using 
equation (4) , for the passive scalar measured in the atmo- 
sphere (inverted triangles) and from a direct numerical 
simulation of three-dimensional homogeneous isotropic 
turbulence 0| (upright triangles). They are virtually 
indistinguishable from each other and from the present 
exponents extracted for thermal convection from fig. 2. 

Noting this good correspondence between z p for ther- 
mal convection and for the passive scalar, and further 
that there is no scaling of unnormalized structure func- 
tions in convection, we suggest that the structure func- 
tions (ATP) in convection should have the general non- 
scaling form 

(AT^) ~ t^ p ip(t) p , (5) 

where the exponents £ p are taken from the passive scalar 
problem. This suggestion is completely consistent with 
the experimentally observed multiscaling properties (2)- 
(4) (fig. 3). One can find f{r) explicitly from the data 
by dividing (AT?) by and taking the p-th root of 
the resulting function. Figure 4 shows the outcome of 
these calculations for different values of p (the calculated 
function <p(r) is normalized by its maximum). In the 
semi-log scales used in the Fig. 4, the ip(r) function has 
an essentially parabolic form given by 

^r) = [l-a(\nr/r h ) 2 } (6) 

where a and r/j are constants. The scatter of the data in 
fig. 4 needs a brief comment. The range of r shown in the 



FIG. 3: Scaling exponents z p , for the normalized structure 
functions S p (r), extracted from fig. 2 (circles). The scaling 
exponents z v calculated using equation (4) for the passive 
scalar in the atmosphere |£j (inverted triangles) and in the 
direct numerical simulations of mixing in three-dimensional 
homogeneous isotropic turbulence [Tof (upright triangles). 

plot covers more than two decades. Although the tails 
(especially for large orders p) do not collapse perfectly, 
most of the data (in the range —2 < lnr < 2.5, i.e. about 
two decades of r) possess a scatter of less than 6%. 

The case of particular interest, corresponding to the 
second order, is shown in the inset to fig. 4; the argu- 
ments due to Corrsin and Obukhov yield a specific 
value for the second order (£2 = 2/3 without intermit- 
tency corrections). The solid parabola (the best fit) cor- 
responds to (6) with a ~ 0.04 and th ~ 0.64s. A slightly 
different choice for (2 to account for the intermittency 
effects shows no essential difference. 

It may be useful to make the following observation. 
In general, for the atmosphere, the structure function 
data are interpreted with respect to spatial separation 
through the use of Taylor's hypothesis. Because of the 
presence of the mean wind in thepresent experiment, 
one may use Taylor's hypothesis y| with equal facility 
and interpret temperature increments over time intervals 
as increments over equivalent spatial distances traversed 
by the mean wind. This hypothesis is neither critical nor 
necessary, though we shall use it later to be in confor- 
mity with standard practice. If the Taylor's hypothesis 
is applied to temperature fluctuations in convection, one 
would replace the time separation r by an equivalent spa- 
tial separation r, and we would have 

(ATP) • [l-o(ln r/mff. (7) 

Local scale invariance characteristic of structure func- 
tions with scaling properties can often be interpreted, if 
only loosely, in terms of a conformal invariance, which 
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FIG. 4: The function tp(r) i n semi-log scales. The solid 
parabola indicates equation (6). 

may be understood as an extension of the Kolmogorov- 
Obukhov similarity hypothesis. There is evidence that 
this conformal invariance is broken for the vel ocity field in 
three-dimensional hydrodynamic turbulence Equa- 
tion (7) suggests the Corrsin-Obukhov scaling possesses 
a symmetry-breaking property (cf also [l2T|). 

To understand the "correction" of the passive scalar 
scaling in (5) it is useful to consider the equations of 
motion in the Boussinesq approximation: 

Pr-\^- + v ■ V)v = -Vp + Tz + Rq- 1/2 V 2 v (8) 



(— +v V)T = Ra-^ 2 V 2 T (9) 

V-v = 0. (10) 

Here v and p are the velocity and pressure fields, z is a 
unit vector in the vertical direction and Pr is the Prandtl 
number. We expect that the anisotropy related to the 
thermal origin of velocity fluctuations is diminished in 
the inertial range of scales with Ra — > oo and so may 
regard, for large Ra, the temperature fluctuations in the 
spirit of a "perturbation" over the passive scalar problem. 

Finally, we wish to address briefly the power spectrum 
of temperature fluctuations. There is no simple way of 
relating the form (7) for the second-order structure func- 
tion to a similarly neat form for the spectrum. However, 
considering that the correction to scaling presumably has 
its origin in a perturbation theory, we will use the same 
form as (7) for the spectrum as well, and replace r by 
(where k is wave number) to write 

E(k) ~ k- 5/3 [l -a{\nk/k h f], (11) 
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FIG. 5: The spectral density of the temperature fluctuation 
(solid curve) against normalized frequency (in log-log coordi- 
nates). The dashed curve corresponds to the logarithmically 
corrected Corrsin-Obukhov spectrum (12) and dotted straight 
lines indicate the Bolgiano-like form E(f) ~ f~ 7 ^ 5 and the 
Corrsin-Obukhov form E(f) ~ ,f~ 5/3 . 

or, for the frequency spectrum, 

B(/)~/- 5/3 [l-a(ln//A) 2 ], (12) 

where fh = t^ 1 and the power "—5/3" is the Corrsin- 
Obukhov scaling for the passive scalar — corresponding 
to ^2 = 2/3 noted above Figure 5 shows the ex- 
perimentally observed spectrum of the temperature fluc- 
tuations (solid curve). The logarithmically corrected 
Corrsin-Obukhov form (12) is shown in this figure as the 
dashed curve. (The fitting range, clearly seen in fig. 5, 
is a bit smaller than for the normalized structure func- 
tions, where a part of "dissipation" range is also cov- 
ered by the fitting.) We also indicate by dotted lines the 
Corrsin-Obukhov form E(f) ~ f~ 5 ^ 3 and the Bolgiano- 
like form E(f) ~ /~ 7 / 5 , mentioned previously in the lit- 
erature in relation to the turbulent thermal convection 
(refs. 0,0- Hi)- However, the Bolgiano-like fit of the 
spectrum has no support in the second-order structure 
function in our case. 

In conclusion, our proposal is that the temperature 
fluctuations in the convection problem can be treated as 
a "perturbation" over the passive scalar problem, and 
that the "correction" is of the non-scaling form shown 
in (7). Empirical evidence in favor of this proposal is 
presented. Developing this suggestion into a full-fledged 
theory is beyond the scope of the present article. 

We thank L. Biven, A. Praskovsky, V. Steinberg, and 
V. Yakhot for useful discussions and help in calculations, 
T. Watanabe and T. Gotoh for providing their paper 
[lfj| before publication, and the referees for valuable com- 
ments. 



4 



[1] K.R. Sreenivasan and R.A. Antonia, Armu. Rev. Fluid 

Mech, 29, 435 (1997). 
[2] L. Kadanoff, Phys. Today 54, 34 (2001). 
[3] E.S.C. Ching, et al, Phys. Rev. E, 67, 016304 (2003). 
[4] D. Biskamp, K. Hallatschek, E. Schwarz, Phys. Rev. E, 

63 045302, (2001). 
[5] J.J. Niemela and K.R. Sreenivasan, J. Fluid Mech. 481, 

355 (2003). 

[6] E.S.C. Ching, Phys. Rev. E, 61, R33 (2000). 
[7] M. Sano, X.Z. Wu and A. Libchaber, Phys. Rev. A 40, 
6421 (1989). 

[8] A.S. Monin and A.M. Yaglom, Statistical Fluid Mechan- 



ics, Vol. 2, (MIT Press, Cambridge 1975). 
[9] F. Schmidt et al., Europhys. Lett. 34, 195 (1996). 
[10] T. Watanabe and T. Cotoh, Statistics of Passive Scalar 

in Homogeneous Turbulence, (submitted). 
[11] V.S. L'vov and I. Procaccia, Europhys. Lett. 29, 291 

(1995). 

[12] D. Muller, Phys. Rev. D, 51, 3855 (1995). 
[13] V. Yakhot, Phys. Rev. Lett. 69, 769 (1992). 
[14] S. Grossmann and V.S. L'vov, Phys. Rev. E., 47, 4161 
(1993). 

[15] S. Ashkenazi and V. Steinberg, Phys. Rev. Lett, 83, 4760 
(1999). 

[16] J.J. Niemela, L. Skrbek, K.R. Sreenivasan and R.J. Don- 
nelly, Nature 404, 837 (2000). 



